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We investigate the influence of critical spin fluctuations on electromagnetic responses in the 
odd-frequency superconductivity. It is shown that the Meissner kernel of the odd-frequency 
superconductivity is strongly reduced by the critical spin fluctuation or the massless spin wave 
mode in the antiferromagnetic phase. These results imply that the superfluid density is reduced, 
and the London penetration depth is lengthened for the odd-frequency pairing. It is also shown 
that the zero field Flude-Ferrell-Larkin-Ovchinnikov pairing is spontaneously realized both for 
even- and odd-frequency in the case of sufficiently strong coupling with low lying spin-modes. 

KEYWORDS: odd-frequency superconductivity, Meissner effect, spin fluctuation, superfluid density, pene- 
tration depth, zero-field FFLO 



1. Introduction 

Odd-frequency superconductivity, whose gap function 
is odd in frequency, is an old but new class of super- 
conductivity. It was introduced originally by Berezin- 
skii for triplet pairing, 1 and by Balatsky and Abrahams 
for singlet pairing. 2 Soon after the Balatsky- Abrahams, 
the possibility of odd-frequency superconductivity (odd- 
u! SC) was discussed in a wide variety of models, e.g., 
the Kondo lattice model, 3-7 the square- 8 and triangular 
lattice 9 Hubbard model, and the t-J model. 10 

The experimental evidence on the odd-cj SC was first 
pointed out by the authors. 11 In Ce-based heavy fermion 
superconductors, such as CeCu2Si2 and CeRhIn 5 , several 
experiments on the nuclear magnetic relaxation rate 12 
and specific heat 15 ' 16 exhibited behaviors of gapless SC 
near the antiferromagnetic (AF) quantum critical point 
(especially in the AF phase). The origin of this gap- 
less SC had long been a mystery. In Ref. 11, the au- 
thors showed that the odd-w SC is realized due to spin- 
fluctuation near the AF quantum critical point and/or in 
the AF state by solving the Bethe-Salpeter equation 17 of 
SC based on the itinerant-localized duality theory. 18,19 
The phase diagram so obtained agrees well with the ex- 
perimental phase diagram, considering the fact that the 
odd-cj SC exhibits the gapless behavior. Accordingly, the 
origin of the gapless SC in Ce-compounds can be under- 
stood as the odd-w SC. 

After this suggestion, detailed investigations have been 
made into the triangular-type lattice, 20-22 and into the 
effective model of U-compounds, 23 both of which are re- 
vealed to possess the odd-w state under some realistic 
conditions. It now appears probable that the odd-w SC 
can be realized in surprisingly familiar situations. The 
concept of the odd-w SC is so general that nowadays it 
has expanded into various fields, e.g., magnetic or den- 
sity wave order, 24 ' 25 superconducting junctions, 26-43 vor- 
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Besides the studies on the possible existence of odd-w 
SC, the thermodynamic stability or the sign of the Meiss- 
ner response had been still remaining unclear. Within a 
naive calculation, the Meissner kernel of the odd-w pair- 
ing seems to have an opposite sing from that of the usual 
(even-w ) pairing, suggesting the thermodynamically un- 
stable superconducting states. Very recently, this puzzle 
was solved by the appropriate treatment of the gap func- 
tion with retardation in the path-integral formalism. 48 ' 49 
According to these works, the odd-cj SC is thermody- 
namically stable, and exhibits the conventional Meissner 
effect. However, these arguments are restricted only to 
the coherent part (without the correlations due to the in- 
coherent part of the quasiparticle) . By contrast, previous 
studies 3-9, 41, 20-23 demonstrated that the odd-cj solutions 
are realized only when the Cooper pairs are mediated 
by a certain strong spin fluctuation, strongly suggesting 
the relevant contribution of the incoherent part. So we 
have to discuss carefully the Meissner effect of the odd-w 
pairing considering the contributions from the incoherent 
part, which are neglected in the previous arguments. 48,49 
In this paper, we consider the incoherent corrections due 
to critical spin fluctuations which are crucial ingredients 
for realizing the odd-w SC. 

In §2, the ordinary arguments for the electromagnetic 
response, where the incoherent part is neglected, are 
briefly reviewed. In §3, we introduce spin fluctuations in 
a general form and investigate the Meissner effect by cal- 
culating the current-current response function with the 
first order correction of the critical spin fluctuations or 
spin wave modes in the AF state. There, we examine the 
singlet pairing in the AF phase, where the realization 
of the odd-w is guaranteed by the previous work. 11 The 
results indicate that the Meissner kernel is reduced due 
to the correction of such spin fluctuations. Next, in §4, 
we study the triplet pairing with ferromagnetic spin fluc- 
tuations in the same manner as in §3. It is shown that 
the Meissner effect is also reduced in triplet pairings. In 
§5, we give another argument of the electromagnetic re- 
sponse by calculating the superconducting density, which 
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is given by the spatially gradient term in the free energy. 
Up to the first order of the critical spin fluctuations, its 
result is equal to that obtained by the current-current 
response function as expected. With this procedure, we 
carry out the calculation of higher order corrections up 
to the third order. Then, we conclude that the Meissner 
kernel is strongly reduced both for even- and odd-w pair- 
ing, even if we take into account the infinite order of cor- 
rections. Our results also suggest that the Flude-Ferrell- 
Larkin-Ovchinnikov (FFLO) pairing is spontaneously re- 
alized without a magnetic field in the case of sufficiently 
strong coupling. Finally, in §6, we discuss implications of 
the present results. 

2. Electromagnetic response of the coherent 
part 

In this section, we reintroduce the Meissner kernel in 
the Ginzburg-Landau (GL) region with the coherent part 
only in a convenient form for later discussions. The rela- 
tion between the current density and the vector potential 
has the following form 50 

j(q) = -Q(q)A(q) 
Ne 2 



rn 



A(q)~ 2 -^Y,P\P- A ^)\ 
m z — 4 

p,n 



x p{p + p(p-) + &{ P+ )& + (p-)\ 
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where p = (p,e n ), p± = p ± q/2, and are the 

normal Green function and the anomalous Green func- 
tion, respectively. N is the number of quasiparticles with 
a charge e and a mass m. The kernel Q xx {q) is given ba- 
sically by the current-current response function, which 
consists of a particle-hole diagram. (Note that the spa- 
tially gradient term in the free energy consists of particle- 
particle diagrams and will be discussed in §5.) These 
Green functions are given as 



A+(p, £n ) 



(2) 



(3) 



e 2 -K 2 + |A M (p,e n )| 2 

both for the even-w (A e ) and odd-w (A Q ). 49 In the GL 
region, i.e., near T c , we can expand them with respect to 
A M . Up to the second order in A M , we have 



Sf (p) ~ G(p) 



|A„(p)|' 



= G(p) + G(-p)lG(p)} 2 \A,(p)\ 2 , (4) 

„ 4 +) (p) 



-G(p)G(-p)AW(p), (5) 



where G is the Green function in the normal state 
G(p) = [ie n — £ p ] . The diagrammatic expressions of 
these expansions are shown in Fig. 1. The Meissner ef- 
fect can be demonstrated in the limit of q — > 0. After 
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Fig. 1. Diagrammatic expressions of the expansion of the normal 
'S and anomalous Green function & '. 



some straightforward manipulations, we have 
QM0) = 2e 2 TY,(^) 2 

p , 71 

x {-2[G(p)} 3 G(-p) + [G(p)G(-p)} 2 } |A M ( £ „)f 



TrV |A ^ £w) 



>0, 



(6) 



which leads to the ordinary (diamagnetic) Meissner ef- 
fect. 50 (Here, we have used the relation N = p^/3tt 2 .) 
This expression, however, does not include the effect of 
the incoherent part, or the spin fluctuations, which is 
indispensable for the realization of the odd-w SC. Only 
from eq. (6), we cannot see how the Meissner effect and 
electromagnetic responses will be modified by the indis- 
pensable corrections of the incoherent part. 

It should be noted here that there is no difference be- 
tween the even- and odd-w pairing in the calculation 
of the Meissner kernel (except for the £ n -dependence of 
A(e„)), since the structures of ^ and & are the same. 
In other words, the classification of even- and odd-cj is 
that with respect to the relative coordinate, p, whereas 
the electromagnetic response is to the center of mass co- 
ordinate, q, so that the theoretical framework of the elec- 
tromagnetic response is basically the same for even- and 
odd-Lu pairings. 

3. Singlet pairings 

3.1 Model 

Our purpose here is to calculate the corrections of 
fluctuations to the Meissner effect. For this purpose, we 
start from a simple and essential model. Judging from 
the previous studies, 3-11, 20,21 the odd-w pairing tends to 
be mediated by spin fluctuations, which originate from 
the local component, i.e., the incoherent part. Therefore, 
we introduce the following low-energy effective action on 
the basis of the itinerant-localized duality theory: 18 ' 19 



S = J2 



1/T 



dT1p p ^ a {T){d T + Cp)'0p : a(T) 
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and g is a coupling con- 



where ^p. a {T) — e 

stant. The spin density is given by 

S(q, T) = 2 i ) P+q/2A T ) (T ^p-q/2A T )' ( 8 ) 
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where a denotes the Pauli matrix. The actual odd-aj pair- 
ing, i.e., the gapless superconductivity, has been reported 
mainly in the coexistence phase of superconductivity and 
antifcrromagnetism. 13-16 For such a situation, the most 
dominant fluctuation is the transverse spin-fluctuation 
due to the spin wave as 

9 2 T N 
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where g is the coupling constant, Tn expresses the energy 
scale of the AF transition temperature, v s is the velocity 
of the spin-wave, and q = q — Q (Q is the magnetic 
order vector). This type of spin fluctuation has already 
been shown to realize the odd-w SC in the AF phase. 11 
When we consider in the doubly folded Brillouin zone, q 
is replaced by q, so that @(q,uj m ) is singular for q ~ 
and uj m ~ 0. We use this transverse spin-fluctuation for 
singlet pairings. 

The corrections of the spin-fluctuation to the kernel 
Qxx are 

AQ xx (q) =2e 2 Tj2@(q')\^(.P)\ 2 

p,q' 

x [-v p v- p &(p + )&(p+ - g')^+(p_)^ t ( P+ ) (a-1 

- v p v^(-p + )&(-p+ + q').J?+( P -)^ lt (p + ) (a-2 
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- u p v_ p Sf(-p_M-p_ + q').^+{p^ lt ( P+ ) (a-4 

+ w p »_ p Sf(p + )Sf (-p+)^+(p_)^ t (-p + + q') (b-1 

+ « p «_ p Sf(p_)sr(-p_)^(-p_ + q'W n {p+) (b-2 

+ WpC-^Sf (p+)Sf (-p- + g , )^(p_)^ t (-p+ + </) 
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Fig. 2. Diagrammatic expressions of the current-current correla- 
tion with correction of fluctuations. 



+ v p v p ^(p + )^{p+ - (p_)^ t (-p_ + </) 

(i-2) 

+ « p » p _,/^+(p_)^ 1 , t (p + )^+(-p + + + 

G-i) 

- « p «_ p ^(p_)^ t (p + )^(-p + + g')^t(P+) (k" 1 ) 

- Up«_ p ^+(p_) t Fi t (p + ) t F+(p_)^ Lt (-p_ + q')] , 

(k-2) 

where v p — p x /m is the velocity of quasiparticles. The 
corresponding diagrams are shown in Fig. 2. In the limit 
of q — > 0, we have 

AQ xx (0) = 2e 2 T5>^(0)|A Al (p)| 2 

3Sf (p)Sf(p)Sf(p)Sf(p) + 35f(-p)^(-p)#-+ (p)^Lt(p) 

+ 2Sf (p)Sf(p) J?+ (p)^t(?) + ^(P)^(p)«^tt(-P)"^t(-P) 
+ 3Sf (p)Sf (p)jr+ (p)^ t (-p) + 3Sf (p)Sf (p)^+ (-p)^t(p) 

- 2Sf(p)Sf(-p)^(p)^ t (p) - 2Sf(p)Sf(-p)^(p)^ t (-p) 

- 2^(p)^(-p)j£+(-p)J%(p) 

+ ^(g , )^(p)^t(p)^(9 / )^tt(-P)"^t(-P) 
+ ^(g')^ t + i(P)^tb)^(9')^t + i(-P)^itb) 
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Here, we have put g' — > since only $>(q, u) m ) with uj m — 
is important in the Ginzburg-Landau region. Up to the 
second order of A(p) (eqs. (4)-(5)), 

AQ xx (0) =2e 2 Tj24^M\ 2 

p 

x [-12G£g1 p + \2G\G 2 _ p - QG 3 p G 3 _ p ] (11) 



2e 2 TVi?o|A Al ( e „) 



,2N F e F 
3m 



{- 12 (S +12 H)- 6 

2e 2 T]>>|A,(e„)| 2 ^^(-6.) 



e 2 iV 



m 



T^^ |A m ( £ „)| 2 (-67t) 



(12) 



where N F = mp F /2-K 2 , e F = p 2 F /2m 1 and G(p)'s are 
abbreviated as G p 's. Note that A(e„) includes the fac- 
tor from the momentum integration of A(p) as A(e„) ~ 
J2 P A(p, e„)$(p, £„), where $(p,e„) is some function, 
and includes the factor from the frequency summa- 
tion (Tj2 q > )- At around T c , the dominant contribution 
of A(e„) comes only from e = 7rT c , so that we can ap- 
proximate as 



% = T V 5 2 7n/« s V ~ ,g 2 T c /T N . 



9 



(13) 



(Here we use u s g c <~ Tn, where g c is a cut off in the 
momentum space.) Finally, the total kernel becomes 



Qxx(o) 



e2N ^v- 

ttT > 



6 T 



|A M ( £ „)| 2 . (14) 



From this result, it is revealed that the Meissner ker- 
nel is reduced by the critical spin fluctuation both for 
the even- and odd-w pairing. When the spin fluctuation 
is moderate, % < (ttT c ) 2 /6 (i.e., g < tt^T c T n /6), the 
Meissner kernel becomes small, but remains positive. On 
the other hand, there is a possibility that the Meissner 
kernel becomes negative when the critical spin fluctua- 
tion is so strong, i.e., 2>q > (ttT c ) 2 /6. However, this rapid 
reduction is relaxed in we take into account the higher 
order corrections as is shown below in §5. 

3.2 Dependences on band structure 

In the derivation of eq. (12), we assume that the den- 
sity of state (DOS) is constant with respect to £ p . How- 
ever, the DOS somewhat depends on £ p in general. Here, 
we see the validity of the above result by expanding the 
DOS up to 0(£ 2 ) as 



E 

p 



v 2 p F(0 = <%N F s F (l + ai£ + a 2 i 2 + ■■■ 

(15) 

The ^-linear term in DOS gives no contribution since 
F(£) is an even function in £ for the relevant quantity, 



cq. (11). For the zeroth order of term (eqs. (6)), the 
contribution from the £ 2 -term gives 



y^yl (-2GpG~ p + G 2 G 2 _ p ) \ ^ 2 ) 
p 



= o, 



(16) 



and for the 0(@q) term, the contribution from the £ 2 
term becomes 



]>> 2 (-12GlGl p + 12G 4 p G 2 _ p - 6G 3 p G 3 _ p ) 



o(? 2 ) 



+ 12x 0-6 



a 2 7r 



0. (17) 



Consequently, the correction to Q xx is exactly zero up to 
£ 2 , i.e., the present result does not depend on the details 
of the band structure. 

4. Triplet pairings 

Next, let us consider the case of triplet pairing. The 
triplet pairing is mediated by the ferromagnetic spin fluc- 
tuation for both the even- and odd-frequency (The de- 
tailed discussions for the relationship between the pair- 
ing symmetry and the spin fluctuation are given in 
Appendix.) We consider the ferromagnetic spin fluctu- 
ation in the form 

g 2 N F K 2 



(18) 



k 2 + q 2 - \u)/r}(q) ' 

where k and kq are the inverse correlation lengths with 
and without magnetic correlations, respectively. n(q) 
is defined as 77(g) = T s fg using a characteristic spin- 
fluctuation temperature T s f. 52 

With the same procedure as for the singlet pairing, we 
obtain the correction to the kernel due to the ferromag- 
netic transverse spin-fluctuation as follows: 

AQ xx (g) = 2e 2 T^(g') 
p,q' 
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- v p v- p+q ,&(p+Wp-Wf i (p+ - q'Wu(p- - q') 

(d-1) 
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{A-2) 

+ V pi V pi <$(p 1+ )&(p 1+ )&( Pl _)&(p 1+ - q') 
+ V p V p &(p+)&(p-)&( P -)<g(p- + q') 

- V p V p &(p-)&(-p+ + q')^+(p + )^ n (p+) 



- v v v. 
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- Up «_ p ^+(p_)^ n (p + )^+(p + - q'W n (p+) (k-i) 

- V p V- p &fHp-)&x(p + )&fi(p-)&x(p- - q') (k-2) 

The corresponding diagrams are the same as that of sin- 
glet pairing except for (b), (c), (g), and (i) type diagrams 
in Fig. 2. (Of course, the spin indices are different from 
the singlet case.) Up to the second order of A(p), 

AQ xx (0) = 2e 2 TY / 4®°\^(p)\ 2 

p 

x [-V2G{pfG{-pf + 6G( P ) 4 G(-p) 2 - 2G(pfG(- P f] 



x{-©-H)-<I)}r 

^^^^ |A,(^)| 2 (-37t) * 



(19) 



Here, % ~ g 2 N ¥ T c nl/ k 2 . The total kernel for triplet 
pairing is then obtained as 



Q xx (0) = — iTV 

m — ' 



l 



|A M (s„)| 2 . (20) 



The Meissner kernel for triplet pairing is also reduced by 
the spin fluctuation both for even- and odd-w , although 
the coefficient of the second term of eq. (20) is smaller 
than that for singlet pairing. 

5. Higher order corrections 

So far, we have shown that the coefficient of AQ XX is 
reduced at least up to the first order of 3l. However, it 
is naturally expected that the signs of the higher order 
corrections are oscillating; when the sign of the first or- 
der is negative, that of the second order will be positive, 
and that of the third order will be negative etc.. Thus, 
we have to check the effect of higher order corrections, 
especially paying attention to their convergence. Unfor- 
tunately, it is too complicated to proceed to the higher 
order corrections with the scheme we took in the pre- 
vious sections. Here, we instead take a strategy where 
we calculate the superfluid density by a different method 
and then determine the Meissner kernel. 

The superfluid density p s is expressed in terms of the 
Meissner kernel as 

N S (T) 



Ps{T) = 



N 



m 



Q, 



(21) 



On the other hand, p s can also be obtained from the spa- 
tially gradient term, |VA(r)| 2 , in the free energy, 53,54 in 
other words, the coefficient of the q 2 -term. This spatially 
gradient term is expressed by the superconducting pair 
susceptibility, K p . p (q), in the form 

^( T ) = ^£IVp)| 2 |^p-p(?)- (22) 
p 

We should have the same results for the Meissner kernel 
through the calculation of K p _ p (q). In this section, we 
first show that the results in the previous sections can 
also be derived from the calculation of K p _ p (q). Then, we 



see the effect of higher order corrections to the Meissner 
effect. 

5.1 Up to 1st order 

The spatial derivative term without the spin fluctua- 
tion is 

d 2 d 2 



dq 2 



4m ^— ' \e n 3 



(23) 



The corrections due to the spin fluctuation, which is the 
same one as in the previous sections, are given by the 
Feynman diagrams in Fig. 3 (a) and (b). Their contribu- 
tions are 

„,2 



AK^ p (q)=Tj2- 



- 2 (6G 5 G_ P 



N ^v- 
- -1— * T /, 

Am Z — / 



[%G p G 2 _ p — 4GpG 3 _ p ) 

3G p G 2 _ p + G p G 3 _ p ) 
6% 



4m 

for singlet pairings, and 



(24) 



AK;_ p (q) =TJ2%^ [-12G 5 G_ P + 6G p G 2 _ p 2G 3 G 3 p ] 
p 

f^Ef^T (25) 
4m ^ \s n 5 

n 

for triplet pairings. Then the total p s (T) up to 0(S> 1 ) 
are 



i 



6^o 



\^{Sn)\\ (26) 



|A M M! 2 , (27) 



for singlet and triplet pairings, respectively. These re- 
sults are exactly equivalent to the results obtained from 
the current-current response function in the previous sec- 
tions, eqs. (14) and (20). 

5.2 2nd and 3rd order 

Now, we have confirmed the equivalence between two 
different procedures for calculating p s (T) up to the first 
order of spin fluctuation, O^ 1 ). With regard to the 
higher order calculations, we can carry out the calcula- 
tions of the spatially gradient term more easily than that 
of the current-current response functions. The results for 
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Fig. 3. Diagrammatic expressions of the pair susceptibility with 
correction of fluctuations. 



singlet pairing are as follows: 



dq 2 



+ 2Gp + G 1 _ p _ + G p+ G 3 _ p _ 



Am 



(28) 



for the second order (Fig. 3(c)-(g)), and 

Mv K ;_ p ( q ) = -Tj2® 2 q 



P,Q 



dq 2 



ZLr p+ Lr_ p _ 

ZU p+ U_ p _ 

rt/~)4 f~l4 

iG p+ G 1 _ p _ 



p+ —p- 



4G 5 , G 



p+ —p- 



N_ v 168% 3 
Am 7 " ^ leJ 9 ' 



(29) 



for the third order (Fig. 3 (h)-(y)). The total singlet pair 
susceptibility up to third order is given by 

O^^E-^-fA + 30 ' " " 



- 150 



18 



(30) 

Rather surprisingly, the coefficients of each term are inte- 
gers (although each term is given as a fractional number, 
e.g., G 5 G 5 _ p = 357r/128|£„| 9 , or G p G e _ p = -7vr/32| £ „| 9 ). 
Moreover, they are almost given as a geometric progres- 
sion, which can be written in the convergent form: 



i 



i 



6(^o/kr, 



l + h{%/\e n \ 2 ) 



|A M ( £ „)| 5 



(31) 

Therefore, p s and so the Meissner kernel is reduced by the 
spin fluctuations, even if we take into account the higher 
order corrections. With this geometric progression form, 
we can go beyond the perturbation theory even though 
we made calculations only up to the third order. The 
results should be valid as far as it converges, so that our 
results (31) is expected to be valid for any magnitude of 
the spin fluctuations, since it is written in a convergent 
form. Then, we also conclude that the Meissner kernel 
can be negative both for even- and odd-w when the effect 
of massless spin fluctuations is sufficiently large. 

Similar estimations of the higher order corrections will 
be given for the triplet pairing just by changing the coef- 
ficients. The result so obtained will be qualitatively the 
same as that for the singlet pairing, namely, the Meissner 
kernel of the odd-w triplet pairing also becomes positive 
with the critical spin fluctuations. 

6. Discussions 

Summarizing the progress hitherto made, we find that 
the Meissner kernel in the GL region is given by 



Qxx — 



2 N ttT 



E 



i - 



6%/\e n \ 



|A M (e„)| 5 



ZU p+ U_ p 

AG p+ G 2 _ p _ 
2G 7 p+ G 1 _ p _ 



ZU p+ U_ p 

4G p+ G_ p 
2G p+ G_ p 

ZLr p+ Lr_ p 



2Gp + G_p 



171 \ £ n\° ^ L 1 + 5^0/ 

(32) 

whose dependence is depicted in Fig. 4. This con- 
clusion will be valid, even if we consider the infinite or- 
der in 3> . The Meissner kernel is positive but strongly 
reduced by the critical spin fluctuations both for even- 
and odd-w pairing. Meanwhile, our results also indicate 
that the Meissner kernel can become negative both for 
even- and odd-cj pairing, when the critical fluctuations 
are sufficiently large. In eq. (31), the dominant contribu- 
tion from the e„-summation comes only from the lower 
e„, since A(e n ) is rapidly suppressed as e n increases 
when the retardation of the interaction is strong. 11,55 
Then, the condition for the sign change is roughly given 
> T 2 (^ 



by 



c) or 9 ^ V^c^n in the case of sin- 



glet pairing in the AF phase. Therefore, we can con- 
clude that with moderate magnitude of spin fluctuations, 
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Fig. 4. (Color online) Relation between the spin fluctuation ®o 
and the Meissner kernel, which corresponds to the super fluid 
density p s or the square of the inverse penetration depth A7 2 . 



< 3) x , the superconductivity with zero center-of- 
mass momentum is stable, but the superfluid density is 
strongly reduced, while with sufficiently large spin fluc- 
tuations, 3>q > S? x , the uniform superconductivity is no 
longer stable. In such a situation of > the super- 
conductivity with Cooper pairs of the finite center-of- 
mass momentum can be stabilized, since the coefficient 
of <7 2 -term is negative, but that of g 4 -term is positive. 
Then, the FFLO pairing can be realized spontaneously 
without the magnetic field. 

In the following discussion, we make some comments 
on the relationship between the present results and the 
previous works. The qualitative tendency of the present 
results are consistent with the standard theory of a su- 
perfluid Fermi liquid, 56 where Al is given by 



X" 2 



4nNe 2 

CO = ?Sr(i 



1 - Y(T) 
1 + %FfY(T) ' 



(33) 



where m* is the effective mass of quasiparticles, Ff is 
the Landau parameter, and Y(T) is the so-called Yosida 
function. It is noted that eq. (33) is also valid for the 
systems without the Galilean invariance. There, the in- 
coherent contribution for the current-current response 
function is not negligible and is incorporated into the 
Fermi liquid parameter Ff , the w-limit vertex , which 
arises from incoherent processes. 56 In the case of cor- 
related electron systems without the translational in- 
variance, such as the heavy fermion systems, eq. (33) 
is modified as the effective mass m* is replaced by the 
dynamic effective mass mj. 57 (md is given as m&jm = 
^}k— fcp uj—Ep ' where S(/c, uj) is the self en- 
ergy.) The spin- fluctuation increases md, so that it re- 
duces A^ 2 . When the system is near the AF critical 
point or the metal-insulator transition, the Drude weight, 
D = (7re 2 iV/TO*)(l + ±Ff), is reduced. 58 " 61 In this case, 
A^ 2 is also reduced, which is consistent with our results. 
However, in the derivation of eq. (33), the possibility of 
the FFLO state is discarded, so that we cannot discuss 
the possibility of FFLO state only from eq. (33). On the 
other hand, in our theory, we directly calculate the co- 
efficient of g 2 -term, which corresponds to the order pa- 
rameter of the FFLO state, so that we can discuss the 



possibility of the FFLO state, resulting in the possible 
zero-field FFLO state. 

Abrahams et al. (and more recently Dahal et al.) in- 
troduced a composite operator which couples a Cooper 
pair to a spin fluctuation and has the same sym- 
metry as the odd-u; pairing. 5 ' 62 They showed that, 
in the composite-operator condensation states, the su- 
perfluid density is reduced in some models of the 
magnon propagators. Although the magnon propaga- 
tor in Ref. 51, T)(q), resembles our spin fluctuation, 
@(q), their contributions to the Meissner kernel are 
different from ours. For example, in the first order of 
T>(q) or @(q), the composite operator theory leads to 
V(q) [-2GpG- p + GpG^p] (eq. (4.3) in Ref. 51), which 
has the form of the first line of eq. (6) supplemented by 
the magnon propagator, whereas C(^ 1 )-term is given 
by 9{q') [-l2G 5 p Gl p + 12G*G 2 _ p - eGlG 3 ^] (the first 
line of eq. (11)) in our theory. Furthermore, in our the- 
ory, the higher order terms are shown to converge having 
an almost geometric progression structure, even though 
the sign of the higher order of @ is fluctuating. How- 
ever, in the composite operator theory, the contribution 
of the higher order term has not been investigated as 
to whether it converges or diverges; it can differ greatly 
from the present results. Nevertheless, both approaches 
share a common aspect in which incoherent degrees of 
freedom beyond the quasiparticle picture (in a sense of 
the Fermi liquid theory) play crucial roles. 

Finally, we comment on the odd-u; in the non-uniform 
system, such as the normal metal / superconductor junc- 
tions. For the Meissner effect of the pairs in the normal 
metal region, we do not need to consider the corrections 
due to the incoherent part, since the spin fluctuations 
are irrelevant there. 

7. Conclusions 

The effects of the critical spin fluctuations, the cru- 
cial ingredient for realizing the odd-u; superconductivity, 
on the electromagnetic response of the odd-u; supercon- 
ductivity have been examined. The Meissner kernel is 
strongly reduced by the critical spin fluctuations both for 
the even- and odd-u; pairing. Thus the superfluid den- 
sity will be smaller and the penetration depth will be 
longer than the case without corrections due to the criti- 
cal spin fluctuations. Moreover, when the effect is strong 
enough, the Meissner kernel becomes negative, suggest- 
ing the zero-field spontaneous FFLO state. This conclu- 
sion has been reached by calculating two different val- 
ues: one is the current-current response function, 50 and 
the other is the spatially gradient term (the coefficient 
of q 2 -term) in the free energy. 53 ' 54 The results, obtained 
by two different methods, completely agree with each 
other at least up to the first order of fluctuation. This is 
valid both for singlet pairing with the antiferromagnetic 
fluctuation in the coexisting phase and for triplet pair- 
ing with the ferromagnetic fluctuation. The validity has 
been also checked for the shape of the density of state 
in the case of singlet pairing, and the present results are 
shown to be invariant, namely, the present conclusion is 
not affected by the detail of the band structure. 



8 J. Phys. Soc. Jpn. 



Author Name 



A(w) 






i1 \ 




Fig. A-l. Illustration of the frequency dependence of pairing in- 
teraction A(u}). 



Fig. A-2. Illustration of the ^-dependent pairing interaction A(oj) 
and the function A(w) cos (ttlo/cuo ) for (a) r\ ~ wo and (b) rj -C cjo- 
The coefficient Ao is given by the integration of A(ui) and Ai by 
that of A(c<j) cos(7rcj/o;o). We can see from these illustrations that 
Ao 3> Ai for the instantaneous interaction and Ao ~ Ai for the 
retarded interaction. 
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Appendix: Possible symmetries of odd-u; pairing 

Here, we discuss the possible symmetries of the odd-u 
pairing on the basis of a general argument. It is found 
that the odd-w pairing becomes competitive with the 
even-w one when the interaction is strongly retarded. 

We assume that the pairing interaction V (q, w) can be 
separated as 

7(q,w) = $(q)A(w), (A-l) 

where the dimensionality of the system is arbitrary. 
Generally, we can assume that A(w) has a maximum at 
lo = 0, and we parameterize its peak width by n as dis- 
played in Fig. A-l. Since A(w>) is an even function, it can 
be expanded as 

. , n An 7TW , 27TOJ , , . 

A(lu) = — + Ai cos h A 2 cos 1 (A-2) 

2 wo wo 

2 r a , * / \ niruj , A . 

A„ = — / duj A(w>) cos , (A-3) 

^0 Jo w o 

with a sufficiently large cut-off ujq. The pairing interac- 
tion is written in the separable form 

A cvcn {e-e') = — + Ai cos — cos — + • • -(A-4) 
2 uj u; 

A odd (e-e') = A : sin — sin— + ••■ , (A-5) 

for the pair scattering (e, — e) — > (e',— e'). The maxi- 
mum eigenvalue is approximately given by Ao/2 for even- 
frequency pairing and Ai for odd-frequency pairing. If 
A(w) is independent of w, that is to say, rj — > 00, Xq > 
and A„ = for all n 7^ 0. When A(w>) is a slowly vary- 
ing function, i.e., r\ ~ wn, Ai becomes finite but very 
small as is seen in Fig. A-2 (a). Therefore, such instanta- 



(a) ferromagnetic (b) antiferromagnetic 




Fig. A-3. Illustration of the momentum dependence of the pairing 
interaction ( I > (q) for (a) ferromagnetic cases and for (b) antifer- 
romagnetic cases. 



neous interaction only mediates the even-w pairing. On 
the other hand, in cases where A(w>) strongly depends 
on w, i.e., r\ <C wn, the function A(w) cos(7rw/ajo) has al- 
most the same form as A(w), so that Ai is nearly equal 
to Ao- Then, T c of the odd-u; pairing can be comparable 
to that of the even-w pairing. In fact, if the momentum 
part of odd-u; odd (q), which is neglected so far, has a 
larger value than that of even-u; cvon (q), it is possible 
that the odd-u; pairing is stabilized prevailing over the 
even-u; pairing. In order to determine which T c is higher 
and what type of gap structure develops below T c , we 
need more data about the model. 

A.l Spin fluctuation model 

Now, we consider the varieties of pairing on the basis 
of the two-dimensional spin fluctuation model. Let us use 

V r B (q,w) = 3$(q)A(o;), (A-6) 

f t (q,w) = -$(q)A(«), (A-7) 

as the model interaction for singlet (V s ) and triplet (14) 
pairing with the momentum dependence of $(q) given 
as 

$(k-k') = </> + 20i [cos(k x - k' x )+cos{k y - k' y )}, 

(A-8) 

where the ferromagnetic spin fluctuations gives 0o > 
4>i > and the antiferromagnetic spin fluctuation 0o > 
-0i > 0. (See Fig. A-3) 
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Table A-l. Possible pairing symmetries and their coefficients of 
the attractive components for ferromagnetic cases, cpi > (upper 
tabic), and antiferromagnetic cases, <f>i < (lower table). The 
frame represents the odd-cj gap. 



We can express $(k — k') in the separable form 68 

$ cvcn (k _ k ,) = ^ + 2 ^( 7fc7fc , + mr]k ,), (A-9) 

$ odd (k-k') = 2^i (sin fc x sin ^ + smk y smk' v ), 

(A10) 

where 7fc = (cosfca; + cosk v )/\/2 and % = (cosfc^ — 
cosk y )/^/2. The functions -jk, Vk, and sinfe^ are irre- 
ducible representations of the symmetry group of the 
square lattice. Then, adding the frequency part A(cj) = 
A.o/2 + X™ A„ cos(n7ro;/cjo), the pairing interactions have 
the forms 

V s c > c (k -k') = 3 [0 O + 2MlkJk> + »)] 

x (A /2 + Ai cose cose'), (ATI) 

V°-°(k - k') = 3 [20i (sin ^ sin k' x + sin k y sin fc^)] 

x (Ai sinesine'), (AT2) 

y t °'°(/c -*;') = - [<t>o + 20i( 7 fe7fe' + %%')] 

x (Ai sinesine'), (AT3) 

V°'°(k — k') = — [20i (sin k x sin k' x + sin fc y sin k' y )] 

x (A /2 + Ai cose cose'), (A-14) 

where e = ne/uo, and V^ e ' e denotes the even-parity and 
even-w singlet pairing interaction, V^ ' denotes the even- 
parity and odd-w triplet pairing interaction, and so on. 

In this representation, the negative coefficient corre- 
sponds to the attractive channel, so that only the ir- 
reducible representation with a negative coefficient can 
become the superconducting gap function. For the fer- 
romagnetic case, i.e., 0o > 4>i > 0, both odd- and even- 
parity triplet pairings are possible, while the singlet pair- 
ing is not. The even-w triplet pairing has the p-wave 
symmetry, and the odd-u; triplet pairing has the s-wave 
symmetry or the d x 2_ y 2-wetve symmetry. The T c for p- 
wave triplet is related to the eigenvalue 0iAo and that 
for s-wave triplet and for d x 2_ y 2-wave triplet pairing are 
related to 7r0oAi and 20iAi, respectively. (See Table A-l 
for "ferro.".) Here, we omit the extend- s-wave pairing 
for the 7 components. Therefore, if Ai is comparable to 
Ao/2, odd-cj s-wave triplet pairing is expected to arise 
since 0o > 4>i in general. The odd-w s-wave triplet pair- 
ing that appears in the triangular lattice 9 ' 20-22 would 
correspond to this case. 



Similarly, for the antiferromagnetic case, — <j>o < <fii < 
0, the anisotropic even-parity singlet pairing, the odd- 
parity singlet pairing, and the isotropic even-parity 
triplet pairing is possible. The T c for the (4,2 ^2 -wave 
singlet is related to the eigenvalue 30iAn and that for 
the p-wave singlet is related to 60iAi. In the case of the 
s-wave triplet, the T c is related to 7r0nAi. (See Table A-l 
for "antiferro.".) Therefore, if Ai is comparable to Ao, 
the p-wave singlet and the d^^-wave singlet pairing 
are nearly degenerate. The p-wave singlet pairing sug- 
gested in Ref. 11 would correspond to this case. 

To summarize, s-wave and d-w&ve triplet odd-w pair- 
ings are possible for the ferromagnetic case, and s-wave 
triplet and p-wave singlet odd-w pairings are possible for 
the antiferromagnetic case. 
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